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I. Introduction 

Ever since Bekenstein and Hawking discovered that the entropy of a black hole is pro- 
portional to its surface area in the 1970s, more efforts, have been devoted to studying the 
statistical origin of the black hole entropy. In 1985, 't Hooft proposed brick-wall method 
and studied the statistical mechanics of a free scalar field in the Schwarzschild black hole 
background [1]. The reason for the introduction of the 'brick wall' method is that the den- 
sity of states approaching the horizon diverges. Thus in order to avoid the divergence in the 
entropy, a cut-off seems has to be introduced, which can be regarded as a 'brick wall'. Math- 
ematically, this method requires the region of non-zero wave function is limited in r + + e 
and L, where r + is the radius of event horizon and e, L are ultraviolet cut-off and infrared 
cut-off respectively, and the relations:^ < r + ,L > r + should be satisfied. Even though it 
seems very arbitrary to introduce a cut-off in the brick wall method, it has been widely 
applied to scalar field and fermion field in various black hole background [2- 11], where it 
showed that the leading term of the entropy for both bosons and fermions is proportional 
to the area of the event horizon. 

Although the brick wall method achieved great success in the study of black hole en- 
tropy, when it is applied to non-spherical black holes such as Kerr, and Kerr-Newman black 
holes, the calculation seems rather complex[12,13]. In particular, the brick wall method 
need one has to adopt small approximation mass and remove the the term which is not 
proportional to area in the integration. Moreover, if one applies the brick wall method to 
Schwarzschild-de Sitter space-time [14]. Different from non-de Sitter space-time, there are 
two event horizons in Schwarzschild-de Sitter space-time. The two event horizons have dif- 
ferent temperatures. Therefore the radiation between them is not in thermal equilibrium. It 
is clear that one should not use the brick wall method which bases on the thermal equilib- 
rium. In other words, the region [r + + e, L] is not in thermal equilibrium. One should notice 
that former work conducted under the brick wall method tells us that the leading term of 
entropy comes from the contribution of the field very close to the horizon. Thus one can 
chose a thin membrane of quantum fields in the vicinity of each event horizon [15]. If the 
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distance from the membrane is 5, e and S should have the same order. That is to say we 
may assume the thickness of the membrane is also e. And then the fields in the membrane 
[r + + £, r + + 2e] can be regarded as in locally thermal equilibrium. In fact, Hawking radiation 
also comes from the vacuum fluctuation in the vicinity of event horizon. Therefore we might 
regard the two event horizons as two independent thermal equilibrium systems and consider 
them respectively. This method is also applied to non-de Sitter space-time. We call this 
method 'membrane method'. The physical picture of this method is very clear. 

On the other hand, the membrane method enables one free from dealing with the an- 
gular modes in the computation as if he/she can separate field equations into two parts: the 
radial equations and the angular equations. That is because just dealing with radial modes 
can give the Bekenstein-Hawking area law. If the worry is about the coefficient, then in the 
membrane method, as in the brick-wall method, one has to adjust the cut-off distance to 
get the usual coefficient, so missing out modes is not going to make matters worse. 

However, an important issue should be pointed out here is the singularities of the met- 
ric of NUT charged spacetime, which is called Misner strings[16]. In order to avoid the 
singularities the time coordinate must be periodic. In the Euclidean section this forces a 
periodicity proportional to the NUT charge that must be matched to the usual periodicity 
requirement following from the elimination of conical singularities in the (r, t) section. Thus 
the NUT charge and the rotation parameter must be analytically continued. Series of papers 
have worked on it [17-22]. But when involving rotation in spacetime, it is not clear that 
the vanishing of the metric function at the horizon yields the same physics its non-Wick ro- 
tated version. Fortunately, R.B.Mann managed to calculate the entropy of Kerr-NUT class 
spacetime by using a boundary counterterm prescription motivated by the AdS-CFT con- 
jecture [23]. We will take the above issues into account in our present paper. Moreover, the 
geometrical properties of several gravitational instantons have been investigated in [24]. As 
it was speculated by some authors that the geometrical properties play an essential role in 
the explanation of intrinsic thermodynamics of black holes and entropy of a black hole can 
be expressed in the following formula: S = Ax/8 [24], in this paper, we will go further to 
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study the entropy of the NUT-Kerr-Newman black holes. Their entropy of spin fields is 
calculated by using the membrane method ( The gravitational field doesn't be taken into 
account because it seems difficult to separate the field equations). The results show that as 
the cut-off is properly chosen, the entropy in the black hole satisfies the Bekenstein-Hawking 
area law. What is more, it seems that the formula S = A\j '8 has its limitations: the author 
only discusses its application to four dimensional sphere-topology black holes. 

Our paper is organized as follows. In the next section, we discuss the singularities of 
the NUT-Kerr-Newman metric and make some mathematic provision. We sperate equations 
in SecIII. Then we will go ahead, in SecIV, with the reduction of the radial part to the one- 
dimensional wave equation through a series of transformations, we calculate the entropy. In 
the last section, we present our conclusions. 

II. NUT-Kerr-Newman metric 
The NUT-Kerr-Newman space-times can be written in Boyer-Lindquist coordinates as 
[25] (here A = 0) 

A - a 2 sin 2 6 (r 2 + a 2 ) asin 2 6 - A (a - (n ~ fflCOse)2 ) 

ds 2 = a 9 Sm ° dt 2 + 2- — ^ a - Utd<t> 

p 2 p 2 

p 2 (r 2 + a 2 ) 2 sin 2 9 - A (a - ia=3^L\ 2 

A p 2 

where p 2 ,p,A, are defined by 

p 2 = p ■ p*, p = r + i (acosO — n) , 

A = r 2 + a 2 - n 2 + Q 2 - 2Mr, (2) 

where M, a, and n,are the mass, angular momentum per unit mass, and the NUT parameter 
(n is also called gravitational magnetic type mass). This metric has a singularity at 9 = 
and 9 = n. The cos 9 term in the metric means that a small loop around the axis does 
not shrink to zero length at 9 = and at 9 = it. This singularity can be regarded as the 
analogue of a Dirac string in electrodynamics, and is called the Misner string. 
The Eucledian NUT-Kerr-Newman instanton is 



ds 2 = 



A + a 2 sin 2 6 
? 



dt 2 + 2 



( r 2 _ fl 2) asm 2Q _ A ( fl _ (n+acosB)^ 



dtd(f) 
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2 f sin 2 6 + A fa - ^±^?Y 
a -^d<p 2 , 



where 
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p = r 



(n + a cos 9) 2 , 



A = r 2 + a 2 - n 2 + Q 2 - 2Mr. 



It was argued by Gibbons and Hawking that there are two kinds of basic gravitational 
instantons based on the dimension of the fixed point set of the continuous isometry group [26]. 
The first kind contains an isolated fixed point. The second type contains a 2- dimensional 
fixed point set is called a bolt. It also has been shown that entropy can be associated with a 
broader and qualitatively different gravitational system, which contains Misner strings[20]. 
Gravitional entropy arises whenever it is not possible to foliate a given spacetime in the 
Euclidean regime by a family of surfaces of constant time r. Such break down in foliation 
can occur if the U(l) has a fixed point set of even co-dimension df < d — 2 (called a nut) [26]. 
If the fixed point set has co-dimension d-2 then the usual relationship between entropy and 
area holds. However, if there exists additional fixed point sets with lower co-dimensionality 
then the relationship between area and entropy is generalized. Following the idea of [23], we 
find the location of the nut is at r = a = r n and the surface gravity is given by 



Here r + and r_ satisfies A = (r — r + )(r — r_) = 0. The regularity in the (r,r) section 
implies that r has period 2ti/k. If we proceeding along this line, we can obtain the entropy 
of scalar field in Euclidean NUT-Kerr- Newman background that it is n{r 2 + — a 2 ), which is 
in agreement with Ref.[23]. 

Here, we turn to calculate the entropy of NUT-Kerr-Newman black holes by the im- 
proved brick wall method in curvature spacetime. The null-vectors of the Newman-Penrose 
formalism [27] we take 



k = 



(r+ — r_) 
2{r 2 + - a 2 ) 
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1 
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We find that the non-vanishing spin-coeffients [27] are: 

A _ „ cos6» 



iasinO 
tt = — = — ; u 

iasinO 



2p 2 p 
1 



2 _*, a = tt-/3*;/3 



1 dA 



T = 



(4) 



4p 2 

Assuming that the azimuthal and time dependence of our fields will be of the form 
%{m<i>-wt) ^ we g n( j faQfc the directional derivatives are 

-A 



5 = m^d, 



D = Pdp = D ,A = 
1 



pV2 



2p 2 ' 

L+,6* = m*' i d li = ^- f =L 



where D n , D+, L n , L+, K, H are defined by 



^ iK ndA 
Dn = dr + ^ + A^> 



-+ 



n dA 

<% + if + ncotO, 
do — H + ncotO, 
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if 
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sinO 



(n — acosOy 
asinO 
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Thus K and H have the relation 



(5) 



(6) 



K — aHsinO = —p 2 uj. 



(7) 



These differential operators satisfy some identities 
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AD n+1 = D n A 
^D+ +1 = D+A. 
(sin9) L n+ i = L n sin9 

(sinO) L^ + i = L^sinO 

(D m \ ( L ^ mas ^ n ^\ ( imasin6\ ( m\ 

III. Spin fields in NUT-Kerr-Newman space-time 
The Maxwell equations in the Newman-Penrose formalism take on the forms[27] 



(8) 
(9) 
(10) 

(11) 
(12) 
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Using Eqs.(4)and (5), then making the transformations 



!>o = $o,0i = 



V2p 



(13) 
(14) 
(15) 
(16) 



and 



02 

we find that Eqs.(13)-(16) become 

A) + ^r $i 
D -— •!• 



(2p*) 



L + 
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(17) 
(18) 
(19) 
(20) 



From Eqs.(18)and(20),$i can be eliminated to give 
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(AD^ + LqLi - 2iup) $o = 



(21) 



Similarly, from Eqs.(19)and(21)there is 



(AD+D + L L+ + 2iujp) $ 2 = 0. 



(22) 



Assuming $ = R+i(j')S+i{9) and $2 = R-i(9)S-i(9) , we can separate the variables of 
Eqs.(22) and (23) to be 



- 2iru) R +1 = XR +1 , 
(AD£Dq + 2iru) R-! = Ai?_i, 
Lq Li + 2u(acos9 — n) S + \ = — XS + i, 
L Lf — 2u(acos9 — n) S-i = —XS-i, 



(23) 
(24) 
(25) 
(26) 



here A is the separation constant. For the Dirac field, the wave equations for a massless dirac 
particles are [28] 

(D + e - p) F l + (5 + 7T - a) F 2 = 0, 
(A'-j/- 7 )F 2 + (<f + /3-<r)F 1 = 0, 
(D + e* - p*) G 2 - (5 + 7T* - a*) G l = 0, 
(A' - p.* - 7*) G ± - (5 + f3* - r*) G 2 = 0, (27) 

where Fi, F 2 ,Gi,G 2 are 4-component spinors,o;, f3, 7, e, /x, 7r, p, r etc are Newman-Penrose 
symbols, and are the complex conjugates of a, (3 etc. 

All the above equations are also separated by using Newman-Penrose formalism. The 
radial equations are given by 



A^D (A5D+R 



2 



\ 2 R 



LtlLi 



Li [LtS_i 
2 V 5 2 



-A 2 S\i, 



-A 2 S_i, 
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(28) 
(29) 
(30) 
(31) 



For scalar field, the separated equations can achieved directly from the Klein-Gordon 
equation 

1 d ( <9$ \ 

' ^r-S =°- ( 32 ) 



yj—g dx^ \ dx u 

The radial equation is 



— A— R(r) + [(r2 + fl2) - ~ am 
dr dr A 



2 



£-A£-i2(r) + lV ' - 7 ""' J fl(r) = A 2 it>(r). (33) 



We can rewrite the above equation with D n , 

AD D+R = X 2 R . (34) 

The radial equations (23),(29)and (34) can be combined into 

AD s D+ 2s _ 1} - 2s(2s - l)iur] R s = \ 2 R S . (35) 

On the other hand, equations (24) and (28) can be written as 

AD+_ S D + 2s(2s - l)iur\ i?_ s = \ 2 R_ S . (36) 

where s is the spin number. It is clear that as s = 0, s = \ and s = l,Eq.(35) corresponds 
to scalar, Dirac and Maxwell field respectively. 

IV. Entropy 

By using the Wentzel-Kramers-Brillouin approximation and substituting equations R s = 
e */M i n t Eqs.(35), the wave numbers are obtained as follows, 



K 2 2s(2s-l) s(2s-l)(s-l) (dA\ 2 X 
k r-~&+ A + A 2 \~a\) "A" ( j 

Considering while s = 0, |, 1 the third term of the equation is zero, we remove it from our 
equation. Therefore we have 

k r = [(r 2 + a 2 )u - am] 2 + 2As(2s - 1) - AA 2 . (38) 

The horizon equation can be written as 



k2 
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A = (r - r+) (r - r_) = 0, (39) 

where r + ,r_, are the radius of the black hole event horizon, the inner horizon respectively. 
We assume that the boson field is in the Hartle-Hawking vacuum state. The spectrum of 
Hawking radiation in NUT-Kerr-Newman space-time can be written as 

N 2 = - (40) 

where N u ,f3 + ,u and Q + are the radiation intensity, the inverse of Hawking temperature, 
the energy of particles, and the angular velocity of event horizon. Moreover, the angular 
velocity, the Hawking temperature of the black hole event horizon are defined by 

+ =lim^ = ~T~T — 2 ' ( 41 ) 

r ^ rH g<M> r% + a 2 

T - - - 1 - ( r +~ r -) (AO) 
+ 2ir P + A7r(r 2 + + a 2 Y { ] 

The (3 + here is in agreement with that of [25] considering A = 0. 

The distribution of particles corresponding to Eq.(40) is given by 

a i = e ^L)±i > ( 43 ) 

where a/ is the number of particles in the / energy level, oji is the degeneracy of / energy 
level; (5 and u>o denote f3 + and Q + . It is clear that Eq.(43) will be meaningless if uj < uo, 
because a h will be negative for bosons in such case. Thus this demands u in Eq.(40) satisfies 
uo — mVL + > for boson fields. However, just as it was pointed out by R. B. Mann in Ref.[23] 
that the Euclidean time variable should be taken to be periodic, so the Lorentzian time 
variable be periodic with period 87m. Therefore solutions to the wave equation must also be 
periodic, and so u in Eq.(33) must be C/(4n), where C is an integer. And also this demands 
in the following calculation the sum of u> should be replaced by the sum of C. Considering 
Eq.(40) we set E = C/(An) —mVt + . Further calculation will show that this does not change 
our results if we take E as the system energy. Thus the wave numbers refer to the horizon 
can be written as 
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k r = — V(r 2 + a 2 ) (E + mtt + - mQ) z + 2As(s - 1) - AA 2 , (44) 

where f2 = r2 " a2 - The free energy at temperature T + of the boson system is given by 

/? + / = -£Jn(l±e-^), (45) 

where + corresponds to fermion field and — corresponds to boson field. 
According to semi-classical quantum theory, there is 

oo 

E — JdEg(E), 

E o 

where g (E) = u' ^jjp is the states density, u/ is the degeneracy of the fields(for scalar field 
and neutrino field, uj' — 1; for Maxwell field, a/ = 2). The states number is 

r (E) = ]T n r (E, X,m) = JdmJ dX^ J k r (E, \)dr. (46) 

The free energy can be calculated as follows 



-P + f = ±JdEg(E)ln(l±e^+ E ), 
o 



„ , oo r + +2e \ ma x A 



J dr J dX jdm±(e?+ E ±l) 1 

r + +e —A ~ 

yV + a 2 ) 2 (E + mtt + - mQ) 2 + 2As(s - 1) - AA 2 , (47) 



where the separation constant A is the angular quantum number which corresponds to I 
in the spherical space-time case. In this case, the range of magnetic quantum number m 
is —A < m < X . X max corresponds to the fact that k r > (while k r = 0, A reach its 
maximum). Considering fermions field and bosons field, the results can be written as 

7 7T 3 (r 2 + + a 2 ) 3 uj' £ 

ff = • "777 • o • -n, (fermions field) (48) 

Jt 180 (3 4 ( r+ _ r _) 2 V 2,KJ J ' K J 

2 7T 3 (r 2 +a 2 )V £ 
fb = ~45 ' P ' ( r+ _ r _) 2 ' ^' ( bosons f idd ) ( 49 ) 
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where e is the ultraviolet regulator, which satisfies < £ <C r + . This manifests that the 
integral over the quantum number m does not diverge, therefore we need not to regularize 
the m integral. On the other hand, the membrane model illustrates that the black hole 
entropy mainly comes from the vicinity of event horizon. Thus we have taken into account 
the following equation in the integration with respect to m, 

Urn n = Q+. (50) 

r— >r+ 

We also used the median theorem in the integration with respect to r, hence e < rj < 2e. 
Eq.(50) is also used in the integration with respect to r. 

We are now ready to obtain the entropy due to arbitrary spin field of the NUT-Kerr- 
Newman black hole from the standard formula 

S = f%. (51, 
As to fermion field, one componential entropy can be written as 

S - 7 ( r + + fl2 ) V e ( 52 ) 

&1/_ 45>' ( r+ -r_) 2 ' ^ ^ 

There are four components of the wave function refer to fermion fields. Therefore the whole 

black hole entropy is given by 

Similarly, the entropy of boson fields can be obtained as 

S - 8 n * ^ + a2 ) V £ (M) 
bb ~45T 3 ' (r + -r_) 2 ' (54) 

Eq.(51) is in same form as that in [18] We choose the cut-off as ^ = 90/5. Here e and 77 in 

Eq.(47) and Eq.(48) are of the same order. Therefore ^ ~ \ = 90(3, then the entropy in 

Eq.(47) and Eq.(48) satisfies the area law 

Sf = I ■ 4tt (rl + a 2 ) J = 7 -A + u', (55) 
S b = l - ■ 4tt (r 2 + a 2 ) J = \a + oo\ (56) 
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where A + is the area of black hole event horizon. The results in Eq.(55) is in agreement 
with our results in Ref.[29] considering A = 0. 

V. Discussion and Conclusion 
We have discussed the issue that arises singularities in NUT-Kerr-Newman spacetime 
and have studied the entropy due to arbitrary spin field in the NUT-Kerr black holes whose 
Euler's characteristic is over two. Our results is in agreement with the results in [23]. Since 
the cut-off was properly chosen, the NUT-Kerr-Newman black hole entropy is identified with 
the Bekenstein-Hawking area law. As the topology of the NUT-Kerr-Newman black hole 
is special, and its Euler's characteristic is greater than two, we can see that the formula 
S = is not applicable to our case. Therefore, it means that this formula in [24] has 
its limitations: the equation does not apply to high dimensional black holes (x (s 2k+1 ) = 0), 
even does not apply to four dimensional black holes. The NUT-Kerr-Newman black hole 
is such an example. Therefore the relations between a black hole topology and its entropy 
need further investigation. In addition, we can see from the results that the electromag- 
netic, Dirac and scalar field entropies of the following black holes: Schwarzschild black hole, 
Reissner-Nordsrtom black hole, Kerr black hole and NUT-Kerr black hole are embodied as 
special cases of the NUT-Kerr-Newman black hole entropy. 
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